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Abstract 

The supersonic diffuser of a Mach 2.68 bifurcated, 
rectangular, mixed-compression inlet was analyzed 
using a three-dimensional (3D) Navier-Stokes flow 
solver. A two-equation turbulence model, and a porous 
bleed model based on unchoked bleed hole discharge 
coefficients were used. Comparisons were made with 
experimental data, inviscid theory, and two-dimensional 
Navier-Stolces analyses. The main objective was to gain 
insight into the inlet fluid dynamics. Examination of 
die computational results along with the experimental 
data suggest that the cowl shock-sidewall boundary layer 
interaction near the leading edge caused a substantial 
separation in the wind tunnel inlet model. As a result, 
the inlet performance may have been compromised by 
increased spillage and higher bleed mass flow 
requirements. The internal flow contained substantial 
waves that were not in the original inviscid design. 3D 
effects were fairly minor for this inlet at on-design 
conditions. Navier-Stokes analysis appears to be an 
useful mol for gaining insight into the inlet fluid 
dynamics. It provides a higher fidelity simulation of 
the flowfield than the original inviscid design, by taking 
into account boundary layers, porous bleed, and their 
interactions with shock waves. 


Nomenclature 

x . axial coordinate, from ramp tip 
y vertical coordinate, from centerline 

he cowl half-height 

p static pressure 

PT total pressure 

Subscripts 

o fieestream 


Introduction 

On supersonic cruise aircraft at Mach numbers 
above 2.0 to 2.2, mixed-compression inlets are used to 
achieve high total pressure recovery and low cowl drag. 
However, mixed compression inlets are susceptible to 
unstarts, where the nomud shock wave is expelled 
forward out of the inlet, total pressure recovery is 
drastically reduced, and severe forces may result on die 
aircraft The angle of attack or Mach number reduction 
the inlet can tolerate before unstarting is a measure of 
the inlet’s operability limit Boundary layer bleed is 
used to condition the shock boundary layer interactions, 
and to help stabilize the normal shock and prevent 
unstarts. Due to die complexity of interactions, (deed 
patterns are typically laid out using empirical 
guidelines. Overall, the design of a supersonic mixed 
compression inlet is a balance between the conflicting 
requirements of internal performance, external drag, 
operability and weight 1,2,3 

Supersonic diffusers for mixed compression inlets 
are generally designed using inviscid analyses. The 
method of characteristics is typically used, because of 
its exact inviscid results and fast turnaround time. The 
actual flowfield differs from die inviscid analysis, due to 
additional viscous effects and interactions. Therefore, a 
number of analysis techniques accounting for these 
effects have been developed, to more accurately predict 
the performance and flowfield. 

In a zonal approach, the flowfield is partitioned into 
distinct areas depending on the dominant physics, and 
each area is analyzed using a separate technique. 
Typically, the inviscid core is computed using the 
method of characteristics, near wall areas are solved with 
a boundary layer code, bleed regions are modeled by 
manipulating the boundary layer profiles, and oblique 
shock-boundary layer interactions are solved using 
control volume analysis. The zonal approach is useful 
for rapidly evaluating a large number of configurations, 
and for determining boundary layer properties in order to 
place bleed regions. However, more complex flows, 
such as corner flows, subsonic regions, separations and 
vortices, are difficult to model properly. Therefore, 3D 
flows, normal shocks, and operability limits are nut 
accurately simulated. 

In reduced or ’parabolized’ Navier Stokes (PNS) 
analysis, the streamwise diffusion term is neglected. 
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For supersonic flow, this allows the solution to be 
marched downstream in space, computing only one 
axial plane at a time. As a result, PNS is 
computationally efficient, and it is well suited for 
modeling supersonic diffusers, including 3D effects. 
However, PNS may not properly simulate inherently 
elliptical flow features such as large subsonic, regions 
and separated flows. Therefore, it may not be suitable 
j for determining operability limits. 

| Although full Navier-Stokes analysis is 

; computationally expensive, its main advantage is 

l generality. Complex flows and interactions such as 

coma flows, glancing shock-boundary layer 
interactions, normal shocks and separated flows can be 
i modeled, at least in principle. Therefore, within the 

f limits of available computational resources, almost all 

aspects of the inlet flow could be simulated, including 
operability limits, terminal shocks, and subsonic 
difFuser flow. General purpose Navier-Stokes solvers 
can be used with only minor modifications, minimising 
the code development and validation effort However, 
many issues remain to be considered, such as truncation 
error, accuracy of turbulence models, and grid 
sensitivity 4 . 

A number of Navier-Stokes solutions of mixed 
compression inlets have been reported in the literature. 
Of particular interest here are 3D Navier-Stokes analyses 
of rectangular supersonic inlets. Rectangular inlets are 
complicated by sidewall and comer effects, that are not 
present in axisymmetric inlets or in foe original 
inviscid design. Shigematsu et al 3 reported a two- 
dimensional (2D) solution of a rectangular inlet in 
supercritical operation, and a 3D solution of the 
supersonic diffuser, both using foe Baldwin-Lomax 
algebraic turbulence model. They obtained qualitative 
agreement with data in the limited comparisons, and 
stated that foe turbulence model would have to be 
investigated to improve foe results. Reddy and Web 6 
showed substantial 3D effects in a Mach 3 rectangular 
inlet, and obtained good qualitative agreement with data, 
using an algebraic turbulence model and a simple 
uniform mass flux bleed boundary condition. It was 
again suggested that a better turbulence model would 
improve simulation of foe comer flows. Fujimoto and 
Niwa 7 re ported a 3D calculation of a rectangular inlet 
including foe terminal shock. Good agreement was 
shown with data, and a vortex was seen along the cowl 
shock-sidewall interaction. Freskos and Fenanhoat* 
produced 2D and 3D solutions of a Mach 1.863 
rectangular inlet including the terminal shock, but dal 
not report any comparisons with data. 

A medium scale model of a Mach 2.68 bifurcated 
rectangular mixed compression inlet with 30% internal 
area contraction was tested in the NASA Lewis 10x10 



foot cross section supersonic wind tunnel (Figs. 1, 2). 9 
Note that foe sidewall leading edge extends slightly 
forwards of foe cowl leading edge. The fieestrea m 
Reynolds number at Mad) 2.68 was 2.5x10® per foot, 
and foe cowl half-height (he) was 10.67 inches. Figure 
3 shows the theoretical shock structure, designed using 
the method of characteristics. The initial ramp shock is 
followed by an isentrppic compression fan, both focused 
on the cowl lip. The initial cowl shock is also 
followed by an isentropic compression fan; the cowl 
shock is canceled at foe ramp shoulder, and foe 
compression fan is canceled on foe curved ramp surface. 
In foe test program, paramedics were performed on foe 
fieestream Mach number, angle of attack, bleed patterns, 
subsonic diffuser vortex generator patterns, and ramp 
positions. Three configurations were selected for 
detailed study, two of which had the ability to self-start 
following an unstart without active controls. One of 
the self starting configurations gave 89 percent total 
pressure recovery and 16 percent distortion at foe 
compressor face station, with a bleed massflow of 6.9 
percent of capture, in critical operation. The spillage of 
about 4.6 percent of capture was much higher than 
expected 

In previous 2D Navier-Stokes analysis of this 
inlet 10 , extensive parametrics were p er fo rmed on 
turbulence models, computational grids and Meed 
models. The solutions were found to be fairly 
insensitive to these parametrics. An algebraic, one- 
equation and two-equation turbulence models were 
examined The two-equation k-e model of Chien 11 was 
rec ommen ded because it is applicable to a wider range 
of flows than foe simpler models, while it is more 
appropriate to use in an engineering environment than 
the more elaborate techniques, such as Reynolds stress 
modeling or large eddy simulation. Four models for 
porous bleed were evaluated. The ‘unchoked bole’ model 
was recommended because it is foe most realistic 
simulation, and because it provides a prediction of bleed 
mass flow rates. Local wall mass flux was determined 
based on discharge coefficients of uneboked bleed holes, 
bleed region porosity, local flow conditions, md 
plenum back pressure. 

In foe present study, foe supersonic diffuser of this 
Mach 2.68 bifurcated rectangular, mixed-compression 
inlet was analyzed using a 3D Navier-Stokes flow 
solver. A two equation turbulence model, and a porous 
bleed model based on unchoked Meed hole discharge 
coefficients were used Comparisons were made with 
experimental data, inviscid theory, and 2D Navier- 
Stokes analysis. The main objective was to gain 
insight into foe inlet fluid dynamics. 
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Numerical Methods 

Hie flow solver used in the present study Is 
NPARC, • Navier-Stokes solver for compressible 
flows. 1 V J The governing equations are the time 
dependent Reynolds averaged Navier-Stokes equations 
with a perfect gas relationship and Fourier’s heat 
conduction law. These equations are discretized in 
conservation law form with respect to general 
curvilinear coordinates and solved with the Beam and 
Wanning approximate factorization algorithm. 
Although the time dependent formulation of the 
governing equations is used, not all portions of tire 
present version of the code are time-accurate, and it is 
intended for steady state simulations. Viscous terms 
were included in all three dimensions, and default 
artificial viscosity settings were used. 

As re com m ended in the previous 2D p arametri c 
study, tire Chien k-e turbulence model and the unchoked 
hole porous bleed model were used. An algorithm was 
implemented for irregularly shaped bleed regions that do 
not follow tiie grid lines, such as the forward sidewall 
bleed. 

The unchoked hole bleed model determines local 
wall mass flux based on discharge coefficients of 
uneboked bleed holes, bleed region porosity, local flow 
conditions, and plenum back pressure. The discharge 
coefficient curves were based on experimental results of 
Syberg and Hickcox 14 . The tangential velocity at foe 
wall was set to zero, and the static pressure and 
temperature gradients normal to the wall were prescribed 
to be zero. Strictly speaking, these boundary conditions 
are not physically correct, but they are consistent with 
the treatment of solid no-slip wails in the flow solver. 
No attempt was made to me d the roughness effect of 
the porous bled surface. 

Across section of the computational grid is shown 
in Figure 4. The grid consists of approximately 
714,000 grid points in 4 blocks (Table 2). All internal 
wall regions were packed to a typical y+ value of 2, as 
calculated from the solution. The external grid was 
fairly coarse. In order to resolve the shock waves more 
crisply, the grid was adapted as much as possible to the 
initial ramp shock and the cowl shock by slanting the 
cross-stream grid lines. The cowl lip was as sharp as 
the grid packing at the wall would allow, about 0.0001 
inch thick. Taking advantage of the two planes of 
symmetry, only one quadrant of the inlet was analyzed. 

All internal walls were assumed to be turbulent, 
and all external walls were treated as slip surfaces. 
Flow conditions were fixed at foe inflow plane and on 
all external planes parallel to the fte es tre em flow. The 
boundary condition at the exit plane was extrapolated. 

Two configurations Were investig at ed. The ’SSI’ 


configuration has higher bleed mass flow rates, and was 
experimentally shown to be self-starting. The ‘NSS’ 
configuration has leas bleed, but was not self-starting. 

The flow solver was run on a Gray Y-MP 
supercomputer. For each configuration, approximately 
10000 iterations were run, with a time step CFL 
number of 0.3. The L2 residual norm bottomed out 
after a reduction of about 3 orders of magnitude. No 
appreciable change was observed in the solution over 
foe last 1000 iterations. 


Results 

Results of foe present 3D analysis are com p ared 
with 2D analysis using the same cross-sectional grid, 
previous 2D analysis in reference 10 using the same 
numerical methods and a finer grid, and experimental 
dataofWassezbaueretal*. 

Mass flow rates and recoveries are shown in Table 
1. Note that in foe 2D analyses, corner and sidewall 
bleeds were neglected, and therefore the total amount of 
bleed was lessened. Bleed mass flow rates were 
p re di cte d fairly well, with the notable exception of the 
forward sidewall, which was substantially 
underpredicted. Spillage was significantly 
underpredicted, although 3D indicates slightly more 
spillage than 2D. Mass flux avenged throat total 
pressure recoveries are comparable for 3D and 2D 
results, indicating that foe sidewall boundary layer 
losses are negligible in the mass averaged sense. 

Bleed patterns and rake locations for the SSI 
configuration are shown in Figure 5a. In the centerline 
surface pressures (Fig. 3b), computational results and 
experimental data both show the cowl shock impinging 
on the ramp slightly in front of the shoulder. Both 3D 
and 2D results Med to predict the high pressure regioo 
on the cowl surface near x/hc * 3.1, but 3D results 
predicted the lower pressures around x/bc*3.3 ~ 3.3. The 
boundary layer profile at the sidewall rake (Fig 3c) was 
not well predicted. For the boundary layer rake just aft 
of the nunp shoulder (Fig 5d) comparisons between 
computational results and data are inconclusive, because 
the accuracy of the experimentally derived total 
pressures may be somewhat questionable. At this rake 
location, the analysis shows the flow to be highly 
nonuniform due to reflected waves. Therefore, it is 
probably not very accurate to use the tingle static 
pressure tap measurement at the base of foe rake to 
convert the rake pitot pressures to total pressures. 

NSS configuration bleed patterns, rake locations, 
ramp and cowl pressures, and boundary layer profiles are 
shown In Figure 6. For this configuration, the surface 
pressures were . qualitatively well predicted. In the 
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boundary layer profiles, the local overshoot above a 
pressure ratio of 1 is believed to be an artifact of the 
grid resolution. It is difficult to properly resolve this 
region with rapid gradient changes near die top of die 
boundary layer, especially in a 3D problem. 

Visualization of the nuierical results for the SS I 
bleed configuration are shown in figure 7. As seen in 
the centerline Mach contours (Fig. 7a), the flow was 
appreciably different from the original inviscid design 
shown in Fig. 3. Sidewall, ramp, and cowl pressure 
contours (Figs. 7b, 7c and 7d) show that bleed regions 
generated substantial waves, a complex interaction 
existed between cowl shock and the sidewall near the 
leading edge, and the ramp isentropic compression was 
steepened near the sidewall. Sidewall surface 
streamlines (Fig. 7e), or simulated ‘oil flows’, reveal 
several interesting features: the ramp isentropic 
compression caused a secondary flow toward die cowl, a 
small separation existed due to the cowl shock - 
sidewall boundary layer interaction near the leading 
edge, and die secondary flow caused the low energy fluid 
to collect near the middle of the sidewall. As seen in 
the detail view of the streamlines near the sidewall 
surface (Fig. 7f), there was indeed a separation vortex 
due to the cowl shock - sidewall boundary layer 
interaction near the leading edge. 


Discussion 

Examination of the computational results along 
with the experimental data suggest that the cowl shock- 
sidewall boundary layer interaction near the leading edge 
caused a substantial separation in the wind tunnel 
model. As a result, the inlet performance may have been 
compromised by increased spillage and higher bleed 
mass flow requirements. This separation was not 
evident from the data alone. The analysis appears to 
underpredict the extent of separation, as suggested by 
several discrepancies with experimental data, and past 
experience with this turbulence model. The spillage 
mass flow was underpredicted (Table 1). A large 
separated region near the cowl-sidewall corner leading 
edge may detach the cowl shock and even force some 
reverse flow out the front of the inlet, increasing die 
spillage mass flow. The forward sidewall Meed mass 
flow rate was significantly un dapre di cted (Table 1) for 
both bleed configurations, indicating that the pressure 
over at least part of the Meed region should have been 
higher, as might be caused by a large sep arat ed region. 
For the SSI configuration, the experimentally observed 
high pressure region at the cowl centerline near x/hc ■ 
3.1 was not predicted (Fig. 5b); however, a large 
separated region would cause a shock wave to propagate 


diagonally downstream toward the centerline, and appear 
as a high pressure region. Furthermore, waves arising 
from the two sidewalls would constructively interfere at 
the centerline, producing an even higher pressure. The 
sidewall boundary layer profile (Fig. 5c) was poorly 
predicted, indicating that the sidewall secondary flows 
such as the separated flow were not accurately 
simulated. Past studies have shown that the Chien k-e 
turbulence model tends to underpredict separation due to 
adverse pressure gradients and shock interactions’* 1 *. A 
turbulence model that more accurately predicts 
separations should improve the accuracy of results. 

Computational results and wind tunnel data both 
show substantial waves in the internal flow that were 
notin the original inviscid design, as seen in the surface 
static pressures (Figs. 5b, 6b) and contour plots (Figs. 
7a-d),. The near sonic internal flow was highly sensitive 
to even small perturbations in the near wall flow. Bleed 
regions were seen to generate waves due to streamline 
turning. At the start of the Meed region, die 
streamlines turn toward the wall, thus producing an 
expansion wave. Over the Meed region, the streamlines 
angle into the wall. At the end of the bleed region, the 
streamlines turn sharply away from the wall as die 
boundary layer thickens rapidly at the start of die solid 
wall, producing a compression wave. Boundary layer 
displacement effects caused the cowl shock to steepen, 
impinge on die ramp forward of the shoulder, and reflect 
off instead of being canceled. 

3D effects were fairly minor for this inlet at on- 
design conditions. Although 3D results were quite 
similar to 2D results, some evidence of 3D flows were 
observed, such as waves g en e rated by sidewall Meeds 
and flow migration on the sidewall. Flow separation 
due to the cowl shock-sidewall boundary layer 
interaction, which was underpredicted, may be 
responsible for additional 3D effects. As the operability 
limits are ap proac h ed, 3D effects will become more 
prominent, as shock-boundary layer interactions become 
more severe, bleed mass flow rates increase and the flow 
approaches sonic conditions. Also, for higher Mach 
number inlets, viscous effects and bleed mass flow rates 
are greater, causing more significant 3D effects. 

Overall, Navier-Stokes analysis ap pea rs to be a 
useful tool for gaining insight into the inlet fluid 
dynamics. Although the analysis does not necessarily 
match the data in all instances, it provides a higher 
fidelity simulation of the flowfiekl than the original 
inviscid design* by accounting for boundary layers, 
porous Meed, and their interactions with shock waves. 
Most of tiie discrepancies with data were traceable to 
particular sources. Therefore, ahen interpreting the 
results, it is important to be aware of the strengths and 
weaknesses of the analysis method. 


Conclusions 


The supersonic diffuser of a Mach 2.68 bifurcated, 
rectangular, mixed-compression inlet was analyzed 
using a three-dimensional (3D) Navier-Stokes flow 
solver. A two-equation turbulence model, and a porous 
bleed model based on unchoked bleed hole discharge 
coefficients were used. Comparisons were made with 
experimental data, inviscid theory, and two-dimensional 
(2D) Navier-Stokes analyses. The major conclusions 
were as follows: 

Examination of the computational results along 
with the experimental data suggest that the cowl shock- 
sidewall boundary layer interaction near die leading edge 
caused a substantial separation in the wind tunnel 
model. As a result, the inlet performance may have been 
compromised by increased spillage and higher Meed 
mass flow requirements. TTiis separation was not 
evident from the data alone. The analysis appears to 
undeqpredict die extent of separation, as suggested by 
several discrepancies with experimental data, and past 
experience with this turbulence model. 

Computational results and wind tunnel data both 
show substantial waves in the internal flow that were 
not in the original inviscid design. These waves were 
generated by interactions of boundary layers, bleed and 
shock waves. 

3D effects were fairly minor for this inlet at on- 
design conditions. 3D results were quite similar to 2D 
results, but some evidence of 3D flows were observed. 
In certain other circumstances, 3D effects may be more 
prominent 

Navier-Stokes analysis appears to be an useful tool 
for gaining insight into the inlet fluid dynamics. 
Although the analysis does not necessarily match the 
data in all instances, it provides a higher fidelity 
simulation of the flowfield than the original inviscid 
design, by accounting for boundary layers, porous bleed, 
and their interactions with shock waves. 
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Fig. 6 NSS configuration results. 
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